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Abstract. The effect of Alfven waves on the thermal instability of the Interstellar 
Medium (ISM) is investigated both analytically and numerically. A stability anal- 
ysis of a finite amplitude circularly polarized Alfven wave propagating parallel to 
an ambient magnetic field in a thermally unstable gas at thermal equilibrium is 
performed, leading to a dispersion relation which depends on 3 parameters, namely 
the square ratio of the sonic and Alfven velocities (/3), the wave amplitude and the 
ratio between the wave temporal period and the cooling time. Depending on the 
values of these 3 parameters, the Alfven waves can stabilize the large-scale pertur- 
bations, destabilize those whose wavelength is a few times the Alfven wavelength 
Xaw or leave the growth rate of the short scales unchanged. In order to inves- 
tigate the non-linear regime, two different numerical experiments are performed 
in a slab geometry. The first one deals with the development of an initial density 
perturbation in a thermally unstable gas in the presence of Alfven waves. The sec- 
ond one addresses the influence of those waves on the thermal transition induced 
by a converging flow. The numerical results confirm the trends inferred from the 
analytic calculations, i.e. the waves prevent the instability at scales larger than 
Xaw and trigger the growth of wavelengths close to Xaw, therefore producing a 
very fragmented cold phase. The second numerical experiments shows that i) the 
magnetic pressure prevents the merging of the CNM fragments therefore main- 
taining the complex structure of the fiow and organizing it into group of clouds 
ii) these groups of CNM clouds have an Alfvenic internal velocity dispersion iii) 
strong density fiuctuations (~ lOpcnm) triggered by magnetic compression occur. 
We note that during this event there is no stiff variation of the longitudinal veloc- 
ity field. This is at variance from the hydrodynamical case for which the clouds 
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are uniform and do not contain significant internal motions except after cloud 
collisions. In this situation a strong density fluctuation occurs accompanied by a 
stationary velocity gradient through the cloud. 

Key words. ISM: instabilities - magnetohydrodynamics - turbulence - Alfven 
waves - thermal instability 

1. Introduction 

It is now well established that the atomic interstellar medium is a thermally bistable 
gas which at a pressure of about 4000 K cm^"^ can be in two different phases namely 
the Warm Neutral Medium (WNM) and the Cold Neutral Medium (CNM) roughly in 
pressure equihbrimn (Kulkarni & Heiles 1987, Field et al. 1969, Wolfire et al. 1995, 2003). 
It is also believed that the atomic gas is strongly magnetized with magnetic intensity 
around 5/iG (Troland & Heiles 1986, Heiles 1987, Heiles & Troland 2005). Although the 
structure of the magnetic field is poorly known, some observational evidences seem to 
indicate that the fluctuating part is at least comparable to the uniform one suggesting 
that magnetized waves may be of great importance for the dynamics of the ISM. Another 
important observational fact is the absence of correlation between the magnetic intensity 
and the density in the interstellar atomic gas. Various works have investigated the effect 
of MHD waves on the dynamics of a polytropic or nearly polytropic (non thermally 
unstable) magnetized gas (Dewar 1970, Goldstein 1978, McKee & Zweibel 1993, Passot 
et al. 1995, Falle & Hartquist 2002, Passot & Vazquez-Semadeni 2003). 

It is thus of great interest to investigate the simultaneous role of magnetic fields 
and of the bistable nature of the flow on the physics of the ISM. Field (1965) considers 
the effect of a transverse magnetic field on the thermal instability and generalizes the 
isobaric criterion. In the context of cooling flows, Loewenstein (1990) studies the thermal 
instability in the presence of a static field. Hennebelle & Pcrault (2000) investigates 
the role of an initially uniform magnetic field, analytically and numerically, when the 
thermal condensation is dynamically triggered. They propose a mechanism based on 
magnetic tension to explain the thermal collapse in a magnetized flow and argue that the 
magnetic intensity in the WNM and in the CNM should not be very different. Piontek 
& Ostriker (2004, 2005) study the magneto-rotational instability in a thermally unstable 
gas. 

Here we investigate the effect of a finite amplitude circularly polarized Alfvcn wave on 
the thermal instability, both analytically and numerically. These Alfven waves, which are 
exact solutions of the MHD equations, are non dissipative and are thus very likely to be 
present in a magnetized gas such as the ISM. Similar waves, although at a smaller scale 
and in a regime affected by dispersive effects, exist in the solar wind upstream of the earth 
Send offprint requests to: P. Hennebelle 
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bow shock (Spangler et al. 1988). In-situ satellite observations have clearly identified 
circularly polarized quasi-monochromatic Alfven wave packets probably generated by 
reflected protons. 

In Sect. 2 we present a stability analysis of a circularly polarized Alfven wave propa- 
gating in a thermally unstable gas at thermal equilibrium and obtain a dispersion relation 
which generalizes the relation obtained by Field (1965). Wo then solve this equation nu- 
merically for various regimes and discuss the consequences for the ISM. In Sect. 3, we 
perform various numerical experiments to confirm the analytic prediction and to inves- 
tigate the non-linear regime. Section 4 concludes the paper. 

2. Analysis 

The equations governing, in the magnetohydrodynamic (MHD) limit, the one- 
dimensional motion of a plasma permeated by a uniform magnetic fleld Bq in a slab 
geometry read 



du du Id / 



dt dx pdx \ Stt 

— +u—-—— (3) 
dt dx 47Tp dx 

(^b d , ^ dv 

_ + (4) 
(f + ' = -(^ - '>^S - + ^M<T)a.T) (5) 

where all fields only depend on the coordinate x and time t. 

The ambient field Bq is assumed to be oriented in the x direction. The velocity field 

V has a component u = Vx along the x coordinate and two transverse components 

combined in the complex number w = Vy + iVz- Similarly, we write the magnetic field 

as 6 = By + iBz, the component Ex = Bq remaining constant. The mass density and 

thermal pressure are denoted by p and P respectively and we assiimc a perfect gas law 

P = — - — where R — — — is the universal gas constant, ks the Boltzmann constant and 
p ruH 

p the mean molecular weight in units of the hydrogen mass mu. Heating and cooling 

processes arc combined in a single net cooling function C The parameter 7 denotes the 

ratio of the specific heats at constant pressure Cp and at constant volume ■ Note that 
R 

Cy = 7 T— . Thermal diffusivity due to neutrals is the dominant one. It is isotropic 

(7-l)Ai 

in spite of the presence of the magnetic field and equal to k{T) = b/iCy'qiT), where 
r}{T) = 5.7 10-^{T/IKY'^ g cm^^ (Lang 1974). 

Throughout most of Sect. 2, thermal diffusivity will be neglected for simplicity, leading 
to assume a vanishingiy small Field length. Its effect will briefly be addressed in Sect. 
2.5 and it is explicitely taken into account in the numerical simulations of Sect. 3. 
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Two main non-dimensional numbers can be defined, the sonic Mach number Mg = 

Vq/cs ratio of a typical velocity Vo with the constant sound speed Cs = \ and the 

Alfvenic Mach number Ma = Vo/ca, where = -Ba;/(47r/>o)^^^ is the Alfven speed of the 
unperturbed system. The plasma beta is here defined hy (3 = 

These equations have exact solutions in the form of circularly polarized plane Alfven 
waves of arbitrary amplitude (Ferraro, 1955). They read bo = B± exp [—ia{kox — coot)] 
with constant density po and temperature To, and zero longitudinal velocity uo = 
0. The transverse wave velocity is related to the magnetic field perturbation by 

1^0 = V±exp[—ia{kox — coot)] = 7==- The polarization of the wave is determined 

by the parameter a, with a = +1 {a = —1) for a right-handed (left-handed) wave. In 
absence of dispersive effect we can, without restriction, take a = —1. Moreover, B± can 
be taken real. The dispersion relation reads Wq = c^kQ. 



2.1. Ambipolar diffusion and wave steepening 

Before considering the idealized configuration of a completely ionized medium with a 
circularly polarized monochromatic Alfven wave train of infinite length, it is important 
to estimate the timescale of two important processes which are not considered in this 
study. The first one is the ambipolar diffusion which takes place in weakly ionized gas 
and which induces energy dissipation due to the friction between the neutrals and the 
ions. The second one is the steepening of the MHD waves that takes place when the wave 
polarization is not perfectly circular or when the wave packet is modulated. 



2.1.1. Timescale of the ambipolar diffusion 

The ambipolar diffusion time is given by (Shu 1992): 

tda = ^2 ' '^^> 

where 7da is the friction coeflacient between ions of density pi and neutrals of density Pn, 
L is the typical scale to be considered and B is the magnetic intensity. 

In the case of a molecular cloud, it has been estimated to 3 x 10^^ cm^ g~^ s~^ (Shu 
1992). In the case of the atomic gas, since the mass of the neutral is twice smaller, 7(ia is 
twice larger, leading to 

lib m 

where (,i is the gas ionization and n„ the particle density of the neutrals. 

In the CNM the ionization is about 4 x 10"'*, therefore at the scale of say 0.1 pc, the 
ambipolar diffusion time is about 2 x 10^^ s or about lOMyr. In the WNM, the density 
is about 1 cm~^, the ionization about 0.1, therefore at a scale of say 1 pc, the ambipolar 
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diffusion time is 5 x 10^'* s. In both cases, ambipolar diffusion operates on a rather long 
characteristic time and can thus safely be neglected. 



2.1.2. Timescale of the wave steepening 

The equation governing the dynamics of parallel propagating Alfven waves in the long 
wavelength, small amplitude limit was derived by Cohen & Kuslrud (1974). A similar 
equation was derived previously by Rogister (1971) from the Vlasov-Maxwell system in 
the context of a collisionless plasma, taking into account kinetic effects such as dispersion 
and Landau damping. The main point is that the evolution of the magnetic field compo- 
nents perpendicular to the ambient field involves a nonlinearity of the form adx (\h\'^b), 
where a is a coefficient depending on /3. Prom this formula, one sees in particular that 
circularly polarized Alfven waves, for which |6| is constant, do not steepen. In the case 
where the wave amplitude is modulated, the characteristic time of steepening is given by 
(using Eq. (40) in Cohen and Kulsrud 1974) 

41-^_B|_ 
3 CA 

Denoting by A the wavelength, by Ic the typical length of the amplitude modulation and 
by Us the characteristic growth rate associated with wave steepening, one gets 



A h} 







(9) 



A{l-P)l,Bl- 

It is to be noted that Cohen-Kulsrud equation is not valid for j3 close to unity, a point 
where sound waves have the same phase speed as Alfven waves and have thus to be 
retained explicitely in the nonlinear dynamics. In general however, we see that the ratio 
LOs/uJo depends quadratically on the wave amplitude and is inversely proportional to the 
normalized scale of amplitude modulation. As an example, for a wave train of amplitude 
ba/Bo = 0.5 and coherence length Zc/A = 5, with /3 = 0.5, one gets LUs/i-^a = 0.075, a 
value about five times smaller than the typical growth rates shown in Fig. 1. This process 
of wave steepening thus appears subdominant with respect to the combined thermal and 
decay instabilities considered in this paper. Moreover, as will be demonstrated by the 
numerical simulations, wave steepening does indeed occur but it does not affect our 
conclusions. 

Another process is at play in real three-dimensional situations, namely turbulent 
cascade, which occurs mainly in transverse directions as a result of nonlinear interaction 
between counter-propagating wave packets. The nonlinear eddy turnover time is expected 
to be smaller than the other processes mentioned above for high enough amplitude and 
situations where Alfven waves propagate in both directions in roughly equal amount. 
When waves propagate in a privileged direction, as in the solar wind, the turbulent cas- 
cade is much less efficient, leaving enough time for the other instabilities to develop. The 
interaction between waves and turbulence in compressible MHD is still a very debated 
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topic (see e.g. Cho & Lazarian 2003). Three-dimensional simulations will have to be 
performed to address the competition between all these phenomena. 

2.2. Derivation of the dispersion relation 

We shall now consider perturbations about the exact solution described above, namely a 
circularly polarized Alfven wave propagating in a uniform gas at thermal equilibrium and 
determine the dispersion relation for normal modes. Such an analysis has been performed 
by Goldstein (1978) for an isothermal gas and by Lou (1996) for a self-gravitating gas. 
Let us write the perturbations as 

p = po + Sp.T =^ To + 6T,v = vo + dv, b = bo + Sb,u = 6u. 

The linearized equations write, denoting by a star the complex conjugate. 



dtSp + pqOxSu = 0, 

podt6u = -dJP - ^d^iboSb* + b*oSb), 
on 

podtSv + dpdtvo + poSud^vo = -^B^d^Sb, (10) 
dtSb + dudxbo + bod^Su - B^d^Sv = 0, 
C,{dtST+{^-l)TodJu) = -((|^)^<5p+ (l^j^JT). 

Taking perturbations in the form 

dp = Spex.p{—iu)t + ikx) + c.c, 

ST = ST exp{—iuit + ikx) + c.c, 

5u = Suexp{—iu}t + ikx) + c.c, (11) 

Sv = Sv exp{—i(LU + LUo)t + i{k + k(:))x) + {Sv )* cxjp{i(LU — LUo)t — i{k — ko)x), 

Sb = Sb exp{—i{oj + oJo)t + i{k + ko)x + {6b )* ex]3{i{(jj — (jJo)t — i{k — ko)x) , 

we obtain the following equations 

— wSp + kpoSu = 0, 

i_ ^ _ ^ ^ + 

-po{uJ + (Jo)Sv + {-(JoSp + kopoSu)Vj_ = -^{k + ko)Sb , 

^_ ^ ^ ^ _ 

-po{uj - aJo)6v + (uJoSp - kQpQSu)Vi_ = -^{k - ko)6b , 

47r 

-LjpoSZ + k-^ipoST + ToSp) + ^k{5b^ + 5b~) = 0, (12) 

-{uj + (Mo)5b^ + B^{k + ko)5u- B:,;{k + ko)5v^ = 0, 

— {uj — ujQ)5b + B±(k — ko)Su — Bx{k — ko)6v = 0, 

iCy{-aj6f + (7 - l)Tok6u) = - (dpCfp + drCSf 
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The fluctuations ST, Sp, and 6b can easily be obtained as a function of du in the 
form 

ujSp = kpoSu, 

^ = n n (^kC4-f ~ l)To + -podpc) , (13) 
°^ ^ ~2 5TT^-Lfc±-7- -± r"=^l 



One then gets the dispersion relation 

^-L \ ( ^\ i ( ^\ k f LO \ LU k 



which can be rewritten 



_ ^fc2^±^ ^ fcM^ " t^" ~!" + ^ , (15) 
UJ + iLUc (w- fc)((5 + fc)2 -4) 

~ c^T-C _ TodrC-podpC I ( dC\ ^ ~ 
where — = 7; , = — = — — = — ( — ) , = t^/wo, k = 

k/ko and A B^/B^. 

Note that in the limit A = 0, Eq. p5|) becomes identical with the dispersion relation 

obtained by Field (1965) whereas if £ = 0, it reduces to the dispersion relation obtained 

by Goldstein (1968). 



2.3. Asymptotic behaviors 

Before numerically solving Eq. p5|) . we consider various asymptotic limits. 
2.3.1. Static magnetic field 

It is possible to recover the dispersion relation in the absence of waves for a situation 
where the ambient field is oblique, making an angle 6 with the x-axis. In this case 

UJ^ Bq 

B'x = Bq cos 6 and one must take the limit cjo ~> with — = ca = —i= cos d and 

fco v47rpo 

A = tan Q. It follows that 



2 7 2 2^ + *^fc K a; sm' 



2 , 



k Cg : — - — o TTi 5 (-^^) 

uj + ujJc oj — k'^v\ cos^ u 
B^ 

where we denote — - — — . 

47rpo 

In the case — ti/2 and close to the threshold (a; ~ ^ kcs), we have 

^ - „2 , „2 ■ (17) 

Since the magnetic field is purely transverse, the magnetic tension vanishes and the 
magnetic pressure adds up to the thermal pressure making the gas more stable. The 
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criterion for thermal instability is simply v\u!c + c^o^h < which was first obtained by 
Field (1965). It shows that a purely transverse magnetic field can suppress the thermal 
instability. This is because, in this geometry, the magnetic pressure is proportional to the 
square of the density and therefore dpPtot can be positive even if SpPthorm is negative. 

2.3.2. Instability thresholds 

This section addresses the neighborhood of the instability, a situation where 3(a^) 0. 

We shall first discuss the case 5R(u;) = 0, corresponding to the so-called condensation 
or entropy mode. In the situation where \LUb\ and \uj^\ are smaller than fc, the characteristic 
cooling time is longer than the Alfven crossing time and the dispersion relation reduces 
to 

CO ~ —I . (18) 

/3(4 -k^) + A^ 

When P is not too small sound waves have time to restore pressure equilibrium while 
the gas cools and one expects that the growth rate will be close to the isobaric one. In 
the opposite case where /3 is very small, the growth rate should be close to the isochoric 
one since Alfven waves are not accompanied by pressure or density perturbations. These 
conclusions are easily recovered form Eq. (|18|l which for very small amplitude (A^//3(4 — 
fc^) <C 1) gives 

^ ^ + ^(4_fc2) - ^) ) ' (19) 

while, in the case where /3(4 — k'^)/A'^ ^ 1 it rewrites 

-»(^+ ^^^^2^ ^ (LJb-^c))- (20) 

In the case where A'^ / f3 ^ 1, more likely to be met in the ISM, the effect of the 
Alfven wave depends both on the sign oiLU^ — uJb and of that of fc — 2. In a typical region 
of the ISM with T~ 1000 K, uJ^ — uJb > and therefore the waves stabilize (destabilize) 
the gas if fc < 2, (respectively fc > 2). 

In the limit fc 3> 1 and/or for yl = 0, one recovers the isobaric growth rate, the 
destabilizing effect of the Alfven waves becoming asymptotically small as fc increases. 

If a}h <C fc ^ 1, we find a growth rate similar to that given by Eq. (|17|1 except for the 
factor 2 that divides v\. This is due to the magnetic tension that tends to unbend the 
magnetic field lines making the stabilization of the magnetic pressure less efficient. 

Two different limits are obtained when fc <^ uj^. When the isochoric criterion is not 
verified and u]/k = a remains finite, a situation where the growth rate vanishes with fc, 
one gets 

(l-a)(«'-/3-) = (l-3a)4^. (21) 

In the absence of waves, the growth rate asymptotically approaches u) = ±(3^ k{LOh/Lu^)'^ 
(Meerson 1996). Equation H21() shows in general, assuming a and A small and thus 
~ PuJh/uJ^ + A'^ /A, that the Alfven waves have a stabilizing effect. 
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When uj^ < 0, the growth rate for k <C \uj^\ remains finite, equal to ui = —ioj^, 
independently of the presence of Alfven waves. 

Thus, in the limit of very short and long wavelengths, Alfven waves do not modify 
the growth rates of the thermal instability obtained by Field (1965). 

For completeness, we now consider the case of adiabatic perturbations, corresponding 
to a wave mode with di.{uj) ^ 0. For this purpose, we restrict ourselves to the case of 
low amplitude waves, i.e A <ti 1 and we set uj = ujr + iui. We therefore have S> w^. 
Moreover, it is assumed that ^ 1 and uJf, <ti 1 so that they can be neglected when 
multiplied by uji. With these assumptions, one finds that to the first order 



The real part, ujr — Cgk simply corresponds to a sonic wave. 

When A = 0, one finds that the conditions for thermal stability of a sonic wave is 
uTc — wi, > 0, sl criterion already obtained by Field (1965). 

When A ^ 0, the stability criterion is modified according to Eq. H22|l . showing that 
the effect of the waves depends in a complex way on /3 and k. 

2.4. Growth rate and physical discussion 

We now numerically solve Eq. (|15|) which can be rewritten as a sixth order polynomial 
containing 4 parameters to be specified, namely (3 and A (characterizing the Alfven wave) 
and LOh and ujc (function of the thermal processes). The value of ujc — lOc/^o represents 
the ratio of the temporal period of the Alfven wave divided by the cooling time and can 
be arbitrarily chosen. However, once uj^ is specified, uJb depends on the thermal function. 
In order to estimate this parameter we use the standard cooling function of the neutral 
atomic ISM (Wolfire et al. 1995, 2003) which is used in Audit & Hennebelle (2005) in 
the thermally unstable regime {n — 3cm^'^ and T ~ 500 K). To obtain the dispersion 
relation we integrate Eq. (|14|l for k/k^ between 0.01 and 100 using logarithmic spacing. 
The roots of the polynomials are obtained using the zroots subroutine (Press et al. 1992). 
Here, we restrict our attention to the unstable branch only. In order to verify our method, 
we have reproduced the dispersion relation for the decay instability of an isothermal gas 
presented in Goldstein (1978). 

Fig. n displays the results of the numerical integration of Eq. p4|l . First panel shows 
the adiabatic case for (3 — 0.9, A — 0.5 (dotted line), (3 — 0.5, A — 0.5 (short dashed 
line), /3 = 0.5, A — 1 (dot-dashed line), f3 — 0.1, A — 0.5 (double dot-dashed line), 
P = 0.1, A = 1 (long dashed line). As expected the circularly polarized Alfven wave is 
unstable (decay instability) in a range of k extending to a few times fco . Both the growth 
rate and the largest unstable value of k, increase with f3^^ and A. 

The second panel shows results for ujc = lOc/wq = 1. The full line corresponds to the 
hydrodynamical case whereas the others correspond to the same values of /3 and A as 




(22) 




Fig. 1. Growthrate as a function of the wave number. First panel is for the adiabatic 
case, second, third and fourth panels display results for various values of Wc/wo- Solid 
line corresponds to the case ^4 = 0, dotted line to (3 = 0.9, A = 0.5, short dashed line to 
(3 = 0.5, A = 0.5, dot-dashed line to (3 = 0.5, A = 1, double dot-dashed line to /3 = 0.1, 
A = 0.5 and long dashed line to (3 = 0.1, A = l. 



in the first panel. Various interesting features can be seen, i) When k ^ the effect of 
the Alfvcn wave is to decrease the growth rate and therefore to stabilize these modes 
with respect to thermal instability. For f3 = 0.5, A = 1 and for f3 = 0.1, the modes whose 
wavenumber is smaller than A" ^ 1 are perfectly stable, ii) The intermediate modes 
(i.e fc ~ fco) are more unstable when f3 is smaller and A is higher. This is due to the 
decay instability that the wave undergoes for these values of k. iii) When k ^ oo, the 
growth rate is independent of f3 and A. These features are in good agreement with the 
asymptotic limit w ~ studied in the previous section. Cases with a larger /3, namely 
1 and 1.5, have also been explored but are not exposed here for conciseness since they 
are not directly relevant for the regions of the ISM we consider in this paper. Although 
the decay instability disappears for /? > 1 (for small enough amplitude), we find no 
qualitative difi'erence with the cases (3 < 1. In particular the intermediate wavelengths 
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are still destabilized by the waves. This is in good agreement with the analytical study 
of the instability threshold (see Eq. (|18|l 'l. 

The third panel shows results for ujc ~ 0.1. In this case the cooling time is 10 times 
larger than the period of the Alfven waves. The dispersion relation is qualitatively similar 
to the previous case. Quantitatively however, the intermediate wavelengths are much 
more unstable than the small wavelengths. This is due to the fact that in this range 
of parameters, the cooling time being larger than the dynamical time of the waves, the 
fastest growing instability is the decay instability. Note that since t^o is 10 times larger in 
this case than in the case displayed in panel 2, the value of lo/loq in the hydrodynamical 
case (full line) is 10 times lower. 

Fourth panel shows results for u)c — 10. In this case, the growth rate of the thermal 
instability is much shorter than the growth rate of the decay instability. Therefore the 
only effect of the Alfven wave is to stabilize the gas with respect to thermal instability. 
This effect increases when (3 decreases and when A increases. 

These results suggest that the presence of non-linear circularly polarized Alfven waves 
in a thermally unstable medium like the neutral interstellar atomic gas can have two main 
effects. If the cooling time is short with respect to the temporal period of the waves, 
then the waves stabilize the gas. Therefore the gas can survive longer in the thermally 
unstable domain leading possibly in the ISM, to a larger fraction of thermally unstable 
gas. Since the short wavelengths are the most unstable the trend is that the CNM is 
very fragmented into several clouds. If the cooling time is larger than the wave period, 
then the decay instability makes the intermediate modes (fc ~ fcg) more unstable. In that 
case, the fraction of thermally unstable gas is not necessarily larger (depending on [3 and 
A) but the CNM should be fragmented in structures having a size of about l/(fco x ^), 
where ^ is the density ratio between the CNM and the WNM. 



2.5. Effect of thermal diffusivity 

Here we briefly consider the effect of thermal diffusivity. When this term is taken into 
account, the dispersion relation H15|l becomes: 

_ ^fc2^ + .(g + g^) ^ fe2^2 S3 + fc^^-3^ + fc ^23) 

UJ + i{u!c + Kck'^) (w - fc)((cj + - 4) 

where Kc = K{To)ujQf3 / {Cy poC"^) and Kb = Kc/j. 

Fig. 121 displays the growth rate obtained from Eq. H23() using the fiducial value of 
k{To) given at the beginning of Sect. 2 and for uJc — 0.3, 1, 3. The different curves are 
associated with the same values of /3 and A as in Fig. For small and intermediate 
values of k the shape and the values of w are very similar to the case of vanishing thermal 
conductivity. As expected however, thermal conduction introduces a cut-off at small scale 
(Field 1965). The value of k for which uj vanishes depends on /? and increases (by a factor 




Fig. 2. Growthrate as a function of the wave number when the effect of the thermal 
diffusivity is taken into account. The curve styles are associated with the same values of 
(3 and A as in Fig. ^ 

2 to 3 for the values considered here) when (3 decreases, confirming the trends inferred 
previously, i.e. the CNM should be more fragmented in the presence of Alfven waves 
leading to smaller CNM structures. 

3. Numerical study 

In order to test the analytic results presented in the previous section and to investigate 
the non-linear regime, numerical simulations are performed in a slab geometry. 

For this purpose we use the ID adaptive mesh refinement (AMR) code presented in 
Hennebelle & Perault (1999, 2000). The AMR technique is very helpful to simultaneously 
resolve the sharp thermal fronts (~ 10^'^ pc) and the larger scale (> 10 pc) involved in 
the problem. The code has been extensively tested and the results in the hydrodynamical 
case have been closely compared with high resolution simulations using a second order 
Godunov scheme. The growth rate for the parametric instability of an Alfven wave in 
an adiabatic gas have been calculated and shown to match, within an accuracy of a few 
percents, the results of the first panel of Fig. ^ 

Two numerical experiments are carried out. First, we consider a situation for which 
the gas is initially thermally unstable and we study the development of the thermal 
instability in the presence of Alfven waves. Although these initial conditions are somewhat 
artificial, they are simple and close to the assumption of the analytic analysis, making 
comparisons easier. Second, we setup more realistic initial conditions corresponding to 
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Fig. 3. Density field for 3 snapshots illustrating the development of the thermal instabil- 
ity in the hydrodynamical case. The density of the perturbation increases until the gas 
reaches thermal equilibrium. 

a converging flow of thermally stable WNM. In this case the ram pressure of the flow 
drives the thermal collapse dynamically. 

3.1. Case of gas initially thermally unstable 

In order to study the effect of the circularly polarized Alfven waves on the development 
of the thermal instability, we start the simulation with thermally unstable gas (ri = 3 
cm""^ and T ~ 500 K) with density fluctuations of amplitude 0.5. The cooling function 
is described in Audit & Hennebelle (2005). The computational domain, which initially 
contains 5000 pixels, has periodic boundary conditions and a length of 3 x 10^^ cm, 
corresponding to 20 wavelengths of the initial Alfven wave. 

Fig. Ul shows three snapshots of the density field in a purely hydrodynamic run (no 
MHD wave is present). The initial perturbation grows as a result of the thermal insta- 
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bility. At time t =2.55 Myr the density is about three times its initial value. At time 
t =3.64 Myr a cloud of CNM having a size of about 0.3 pc, has formed. 

Fig. 21 displays the density and the y-component of the magnetic field for 3 snapshots 
showing the development of the thermal instability in the presence of Alfven waves of 
amplitude B± = 1 fiG and for B^c = 5 fiG. In this situation the value of lOc/^o is about 
0.03, A = 0.2 and (3 = 0.17. The peak density at time t = 1.82 Myr is about 4.9 whereas 
the initial peak density is 4.5. In the hydrodynamic case the peak density at the same 
time is about 6.2 showing that the waves have significantly slowed down the growth of 
the perturbation, by a factor of about (6.2 — 4.5) / (4.9 — 4.5) ~ 4. Note that interestingly 
enough, this factor is significantly larger than what is predicted by the linear theory 
which predicts a difference of about 10% (in the case of a perturbation having an initial 
amplitude of 10~^ we verified that the growth agrees with the linear theory). After the 
central density has increased by a factor of about 2 (panel 2), the waves drastically change 
the structure of the gas and create significant density contrasts (~10%). The resulting 
cloud (panel 3) contains density fluctuations of about ~ 50% its maximum value, and 
is therefore very different from the uniform cloud formed in the hydrodynamical case. 
These large fluctuations are due to magnetic pressure variations. It is worth noting that 
due to the contraction, the waves inside the growing perturbation (panel 2 and 3) have 
a larger amplitude and a shorter wavelength than the waves in the surrounding medium. 
According to the analysis of the preceeding sections, this effect tends to increase the 
influence of the waves. 

Fig.Elshows results for B±^ = 2.5 jiG and for B^ — 5 fiG. The waves strongly influence 
the gas evolution. The initial perturbation is totally stabilized by the waves and does not 
develop (second panel). On the contrary, the waves trigger the formation of structures 
having a wave number fc ~ fcp (flrst and second panels) . These structures keep condensing 
(panel 3) and finally about 12 CNM clouds form. Their size is about 0.03 pc and therefore 
about 10 times smaller than the size of the cloud that forms in the hydrodynamical case. It 
is worth noting that the formation of these clouds is significantly (50%-100%) faster than 
in the two previous cases. Therefore in this range of parameters, the waves destabilize 
the gas and accelerate the formation of CNM structures. 

Fig.Elshows results for B± = 2.5 fiG and for B^ = 2.5 fiG. In this situation the value 
of LOc/^o is about 0.06, A = 1 and (3 = 0.34. As can be seen on first panel, the initial 
perturbation does not grow (the peak density is smaller at time t = 1.82 Myr than the 
initial value). Unlike the previous case, the Alfven waves do not efficiently trigger the 
formation of structures at fc ~ fcg (see first panel). However since the waves are unable 
to stabilize small wavelengths, the initial perturbation breaks down in several structures 
which finally develop (panel 2). This nevertheless occurs after a significant delay showing 
that in this situation, the gas spent a longer amount of time in the thermally unstable 
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Fig. 7. Spatial zoom showing the thermal condensations which have been induced by 
the large scale converging flow in the hydrodynamical case. Two snapshots are displayed. 

state. At time t = 5.41 Myr about 6 small structures of size ~ 0.03 pc and one larger 
structure of size ~ 0.3 pc have formed. 

3.2. Case of thermal condensation dynamically induced 

We now consider a converging flow of WNM (Hennebelle & Perault 1999) in a simulation 
box of length 150 pc. The peak velocities are 2.2 Cwnm and -2 Cwnm and the peak to 
peak distance is about 60 pc. Two simulations are performed. The first one is purely 
hydrodynamical, while the second one starts with an Alfven waves of amplitude 2.5 /iG 
with 100 spatial periods in the integration domain. The total magnetic intensity is 5 /xG 
and the longitudinal one — 4/iG. An initial resolution of 25,000 pixels is used in order 
to ensure an accurate description of the Alfven wave (it corresponds to 250 pixels per 
period). 

Figd shows the hydrodynamical result. At time t = 6.09 Myr, two clouds having 
a size of about 0.1 pc have formed. They present weak density gradients and have no 
significant internal velocity. The two clouds have a relative velocity of about 0.15 km/s 
and undergo a collision at time t — 6.3 Myr. A shock-compressed layer forms with a 
density of about ~ 700 cm^'^ and a length ~ 0.02 pc. During the time of the collision, 
the structure presents a stiff velocity gradient. 

FiglH displays four snapshots of the longitudinal velocity field, the density and the 
y-component of the magnetic field in the MHD case. The results displayed in the first 
panel confirm the trends observed in the numerical experiments of Sect. 3.1. The gas 
fragments into few small CNM structures having a physical length as small as few 0.01 
pc (note that due to the AMR scheme, these structures are well described). This situation 
is very different from the hydrodynamical case for which the structures arc much larger 
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and uniform. A comparison between the 2 times displayed in first panel of Fig. |H| reveals 
that all the structures do not form at the same time. The structure of the magnetic 
field deserves attention. It varies very rapidly in the new born condensations {x = 242 
cm and x = 243 cm), therefore compressing them. It varies less stiffly near older struc- 
tures {x ~ 244 cm) since the field lines had time to unbend. Between the clouds, the 
magnetic field is much more uniform. This is a consequence of the fact that the Alfven 
speed is about 10 times larger in the WNM than in the CNM. It is also clear that the 
intercloud magnetic pressure plays an important role in preventing the merging of the 
CNM structures therefore maintaining the complexity of the flow and organizing it into 
groups of structures rather than into a single cloud. In these groups of structures the 
longitudinal velocity dispersion is not negligible, unlike in the hydrodynamical case. It is 
about 0.1-0.2 km/s, i.e. few times the sound speed or the Alfven speed of the CNM. This 
velocity dispersion is due to the transfer of magnetic energy into longitudinal motions 
because of the magnetic pressure fluctuations. We stress the fact that observing such 
group of structures with a low spatial resolution (~0.2 pc) may lead to a rather different 
picture, namely a broad, uniform and turbulent (having a Mach number M~ 1 — 2) CNM 
structure. 

The second panel of FiglHl shows two later snapshots of the same numerical experi- 
ment. The group of structures which is seen in the first panel is now located at x ~ 251.5 
pc. Large fluctuations (~ 100%) of density and magnetic field are still present as well 
as a longitudinal velocity dispersion of about 0.1 — 0.2 km/s. Another smaller group of 
structures {x ~ 253.5 pc) has formed. At time t — 9.11 Myr (dotted line), the first cloud 
undergoes a large density fluctuation (rimax — 600 cm"'^) at a scale of about 0.01 pc. At 
the same time the cloud is compressed, because of the (magnetic) interaction with the 
other cloud, so that its length is divided by a factor ~ 2. An interpretation based on 
the previous analytic results is that the amplitude and wavcnumbcr of the Alfven wave 
in the cloud increasing as a result of the compression, the decay instability is triggered 
leading to a larger density fluctuations. Interestingly, in contrast with the stiff and large 
fluctuations undergone by By, the longitudinal velocity field remains relatively smooth. 

4. Discussion and conclusion 

The analysis presented in Sect|21as well as the numerical experiments discussed in the 
previous section show that even modest amplitude Alfven waves may have a strong impact 
on the structure of the multiphase ISM. This study focused on circularly polarized parallel 
propagating Alfven waves, mainly to allow analytic calculations. These waves, which are 
exact solutions of the MHD equations, are very weakly dissipative and are therefore very 
likely to be present in the ISM. 
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Their effects depend on (3, uTc and A, respectively the square ratio of the sound to 
the Alfven speeds, the ratio between the wave temporal period and the cooling time and 
the wave amplitude. Depending on the values of these parameters, these waves may: i) 
stabilize the wavelengths larger than that of the Alfven wave, Aaw, that would otherwise 
be thermally unstable, therefore enhancing the fraction of thermally unstable gas in the 
ISM , ii) destabilize the wavelengths comparable to Aaw and thus fragment the CNM into 
several spatially correlated small clouds , iii) induce strong density fluctuations within 
preexisting CNM structures (up to 10 times the mean density), iv) maintain an Alfvenic 
velocity dispersion within the CNM structures by pumping their energy into longitudinal 
motions. Finally magnetic pressure tends to prevent the merging between CNM clouds. 

These effects are not observed in one-dimensional purely hydro dynamical simulations 
leading to larger and almost uniform structures with weak internal motions. In two di- 
mensions (Koyama & Inutsuka 2002, Audit & Hennebelle 2005) the situation is more 
complex due to the role of turbulence, with coexistence of thermally unstable gas and 
small-scale structures, but the latter are still relatively uniform and do not contain sig- 
nificant velocity dispersion. According to the present analysis, MHD waves may enhance 
the thermal fragmentation found in these simulations, maintain an Alfvenic velocity dis- 
persion and generate large density fluctuations within the CNM structures. 

Recent observational progress has revealed the presence of interesting features and of a 
large quantity of thermally unstable gas (e.g. Heiles 2001, Heilcs & Troland 2003, Miville- 
Deschencs et al. 2003). Heiles (1997) summarizes and discusses observations of tiny small- 
scale structures and more recently Braun & Kanckar (2004, 2005) and Staminorovic & 
Heiles (2005) report the detection of very low column density CNM clouds. 

Qualitatively, the present study as well as the 2D simulations of Koyama & Inutsuka 
(2002) and Audit & Hennebelle (2005) show similar features including thermally unstable 
gas (also observed by Gazol et al. 2001 in simulations at larger scale), low column density 
structures (down to and even smaller than 10^* cm~^) and large density fluctuations 
at small scale (up to 10^ cm~^). Small-scale structures and large density fluctuations 
appear to be a natural outcome of the two-phase nature of the flow. In the absence of 
turbulence or magnetic field, CNM structures have a typical length of about 0.1 pc (or 
a column density of ~ 10^® cm~^) (Audit & Hennebelle 2005). In the presence of waves 
or turbulence, a growing structure can however fragment into smaller pieces. Moreover, 
since both the sound and the Alfven speed change by a factor ~ 10 between the WNM 
and the CNM, with a transition occurring on a short distance of the order of a Field 
length, supersonic or superalfvenic CNM motions are expected to be more frequent. 

No systematic and quantitative study attempting to closely compare simulations and 
observations has been carried out yet. This is clearly a major challenge for the future. 
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Fig. 8. Same as Fig.[2|in the presence of ch'cularly polarized Alfven waves (see text for 
detail). Four snapshots are displayed. 
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